HIGH-FREQUENCY AVERAGING IN 
SEMI-CLASSICAL HARTREE-TYPE EQUATIONS 
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JOHANNES GIANNOULIS, ALEXANDER MIELKE, AND CHRISTOF SPARBER 

Abstract. We investigate the asymptotic behavior of solutions to semi-clas- 
sical Schrodinger equations with nonlinearities of Hartree type. For a weakly 
nonlinear scaling, we show the validity of an asymptotic superposition principle 
for slowly modulated highly oscillatory pulses. The result is based on a high- 
frequency averaging effect due to the nonlocal nature of the Hartree potential, 
which inhibits the creation of new resonant waves. In the proof we make use 
of the framework of Wiener algebras. 



1. Introduction and main result 

In this work we are interested in the asymptotic behavior for < e < 1 of the 
following nonlinear Schrodinger equation 



where (t,x) G M. + x M d , d G N. This model describes the time-evolution of a 
complex- valued field u e (i, x), under the influence of a Hartree-type nonlinearity (cf. 
[17] for a broader introduction) Above, K = K(x) denotes a given real-valued 
interaction kernel, to be specified in detail below, and " * " denotes the convolution 
w.r.t. x G M. d . The scaling of (1.1) for small e > corresponds to the semi- 
classical regime, i.e. the regime of high-frequency solutions u e (t, x), which can be 
approximated via geometric optics. 

The asymptotic behavior of (1.1) as e — > 0+ has been studied by several authors, 
mainly for the case K{x) = ±jij : For example in [14, 16, 18], Wigner measure tech- 
niques are invoked, which however require mixed states and thus can not be applied 
to our situation. In the one-dimensional case, this constraint can be overcome [19], 
but uniqueness of the limiting Wigner measure for t > is still open. Turning 
to multi-scale WKB expansions, which are typically valid for short times only, a 
variety of asymptotic results can be found in [1, 7, 10, 15], provided that is given 
as a single-phase WKB initial data, i.c 



with given (e-independent) phase f(x) G C and amplitude a E (x) G R, such that, 
asymptotically, a £ ~ a n + sai + e 2 a2 + . . . . 
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In the present work we are interested in generalizing these studies to the case of 
(a superposition of) several WKB waves. Due to the expected nonlinear interaction 
between high-frequency waves (i.e. the appearance of resonances), this problem is 
notoriously difficult, even on a formal level. We shall therefore simplify the situation 
considerably by turning our attention to the case of (asymptotically) small initial 
data corresponding to modulated plane-waves. More precisely, we consider Ug(rr) = 
e a ' 2 u e ] (x), where a ^ 1 and Uq(x) is given by a superposition of e-oscillatory plane- 
waves, i.e. 

J 

(1.2) ui{ x ) = Y. a j( x y k3 ' x/E > k J e Rd ' 

3=1 

with amplitudes aj(x) £ C, j £ {1, . . . , J} C N. Here, and in the following, we shall 
assume for the sake of simplicity that the initial amplitudes aj do not depend on e, 
since we shall only be interested in the leading order asymptotics. A generalization 
to amplitudes admitting an asymptotic expansion in e is then straightforward. 
Without restriction of generality we shall also assume that 

kj ^ k e for all j ^ £ £ T := {1, . . . , J} C N, 

where J E N U {oo} and {1, . . . , oo} means simply N. 

By rescaling u e = £~ Q / 2 u e , we can rewrite the considered model into 

(1.3) ied t u E = -yAu e + e a (K * \u E \ 2 ) u e , u E \ t=Q = u e , 

where the initial data is now of order one but the equation displays an (asymp- 
totically) small nonlinearity. Wc shall from now on take on this point of view 
since it looks more natural from the physical point of view. In addition, it is well 
known (see e.g. [3, 4]) that the choice a = 1 for (1.3) is critical as far as semi- 
classical behavior is concerned (see Section 2 for more details). We shall therefore 
pay particular attention to this case. We remark that the same asymptotic scaling 
has been used in [5, 6] for Schrodingcr equations with (gauge invariant) power-law 
nonlincaritics cx |u| 2cr u, a G N. In particular, in [5] the problem of high-frequency 
wave interaction is exhaustively studied in the case a = 1 for which a geometric 
description of all possible nonlinear resonances is given. 

However, for the considered case of Hartree type nonlinearities, the situation is 
very different, due to the nonlocal nature of K * \u £ \ 2 . Since we expect the solution 
u s (t,x) to be given asymptotically as a superposition of highly oscillatory waves, 
we clearly can not regard the Hartree term to be slowly varying (as in the case of 
a single wave). The notion of a resonance though, seems to be not clearly defined 
in this case. A more sophisticated analysis of the oscillatory structure of K * |u e | 2 
is needed therefore. As wc shall see, the nonlocal nature of the Hartree potential 
yields a kind of averaging effect. In particular no new highly oscillatory phases are 
created in leading order (via resonances) , in sharp contrast to the situation for local 
nonlincaritics. 

In order to be more precise, we need to introduce some notation: The Fourier 
transform of / £ L 1 D L 2 (R d ) will be denoted as 

(mo^m = j^lme-^dx. 
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Our analytical approach will be based on the use of the Wiener Algebra (see Section 
3 for more details). Within this framework it turns out that the natural space for 
the amplitudes a = (a 3 )j 6 r is given as follows. 

Definition 1.1. 

A(R d ) := {a = ( aj ) jer : (flj) j& G * 1 (T;L 1 (R' , ))} J 
r = {1, . . . , J} C N, equipped with the norm 

j 

ll a IU(K d ) = X! H a illi 1 (R d )- 

3=1 

We are now in the position to state the main theorem of this work. 

Theorem 1.2. Ford ^ 1, consider the Cauchy problem (1.3) with a ^ 1, subject to 
initial data Uq of the form (1.2), where the initial amplitudes a,j G L 2 (M. d ) n A(M. d ) 
are such that (dP.aj)j£r G ^4(R d ) /or all \p\ ^ 2. In addition assume 

(1.4) lAI^ 1 := inf{|fc - fc m | : £,m g T, l± m} > 0, 

and let the interaction kernel K satisfy (1 + |£|)if(£) G L oa (M. d ). 

Then, for all T > there exists C, £o > 0, such that for any e G (0, So), the exact 
solution to (1.3) can be approximated by 

(15) lltt'-ti* II <Ce^ with 8 -i 1 ,/a£{1}U[2 ' 00) ' 

77ere f/ie approximate solution u| pp «s given by 

j 

(1.6) <p P (t,^) = X)oj(rc - a-^^V^^-^ 12 /^, 



where Sj(t,x) G R is defined in (2.10) if a ~ \ and Sj(t,x) = if a > 1, respec- 
tively. 

In (1.6), the total number of (highly oscillatory) phases J G N U oo, is the same 
as for the initial data (1.2). Hence, no new phases are created by the nonlincarity. 
Nonlinear effects only show up in leading order via self-modulation of the amplitudes 
(provided a = 1). The condition (1.4) can be seen as a small divisor assumption 
required in the case of infinitely many phases, though, not of the same type as the 
one used in [5]. Obviously, (1.4) is always satisfied if J < +oo. 

Note that if a G N, the error of our approximation is at least of the order 0(e) 
and thus smaller than the size of the original (small amplitude) solution. However, 
we see from (1.5) that if a G (1,2) the error becomes worse as a — > 1+. This 
problem can in principle be overcome if one allows the leading order amplitudes to 
be mildly dependent on e themselves, cf. Remark 2.2 below. 

Finally, one should note that under the general assumptions of this work, we can 
not infer global well poscdness of equation (1.3) in the Wiener space. The usual 
arguments for proving global existence (see e.g. [17]) involve the conservation of 
the L 2 norm, which also holds in our case. However, this is not sufficient to control 
the nonlincarity (K * \u\ 2 )u in the Wiener space. Nevertheless, the above theorem 
shows that for initial data Uq in the form (1.2), the solutions cannot blow up too 
fast: If T E > denotes the blow-up time, then T e — > +oo as e — > 0+. 
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Remark 1.3. A particular example for K, satisfying the assumptions is given 
by the one considered in [2, equation (15)] and more importantly by the so-called 
Yukawa potential 

e -A|x| 

K(x)=± — — — , iel 3 ,A>0, 
\x\ 

for which the corresponding Fourier transform is found to be 

1 



m) = ± 



A 2 + iei 2 ' 



Clearly, (1 + \£\)K € L°°(R 3 ) in this case. Note that in the limit A^Owe obtain 
the Fourier transform of the Coulomb potential K(x) = ±yiy, which, however, is 
too singular, to directly apply our theorem. It remains an interesting open problem 
to establish the same result for the Coulomb case in d = 3. 

The paper is organized as follows: In Section 2, we formally derive the approxi- 
mate solution and draw some further conclusions from it. In Section 3 we set up the 
Wiener framework for the exact and the approximate solution. Finally, we prove 
the required estimates on the remainder of the approximation and consequently 
state the proof of Theorem 1.2 in Section 4. 



2. Derivation of the approximate solution 
We seek an approximation of solutions to (1.3) in the form 

./ 

(2.1) ui PP (t,x) = Y,Mt^y pAt ' x}/e - 

3=1 

Assuming for the moment sufficient smoothness for Aj and </>j, we can plug this 
ansatz into (1.3), which yields 

2 2 

(2.2) fecK pp + £ -Aul pp - e a (K * |< pp | 2 ) < pp = ]T e n X n + e a (Y + Y R ), 

n=0 

where we denote 

J , 

3=1 

J , 

Xx = iJ2 ^ lle (p t Aj + VAj ■ V4> 3 + -AjA^ , 



3=1 

and also 



1 3 

(2.3) X 2 = -Y,e^ /£ ^A 3 



3 = 1 
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These terms are the same as in the linear case K = 0. Due to the presence of the 
Hartree type nonlinearity, we also obtain 

(2.4) Y=-Y,^ l/e (K^\M 2 )A ] , 

j=l l=\ 
,7 J 

(2.5) Y R = ~Y. e ^' e { K * E (Mme'^'-^'^A, 

j = l t,m=l 

Obviously, Yr carries high-frequency oscillations, which are not captured by our 
ansatz (2.1). Thus we need to develop a more careful analysis in the following, 
which shows that Yr is of higher order. 

Ignoring this problem for the moment, we consequently aim to eliminate equal 
powers of e. Hence, in leading order, we set Xq = 0, which is equivalent to the 
Hamilton- J acobi equation 

(2.6) 9t j + I|v^| 2 = O, 4>j\ t=0 = k j -x. 



Solutions to (2.6) determine the characteristic high-frequency oscillations present 
in u| pp . In our case, they are easily found to be 

(2.7) <t >j {t,x) = k j -x- t -\k j \ 2 . 

These phases obviously solve (2.6) for all (t,x) elx R d , i.e. no caustics appear 
in our study. 

In the next step we set X\ = if a > 1 and Xi+Y = 0ifa; = l (note that 
we do not include Yr here). Comparing the prefactors of the terms multiplied by 
e'^/ e 5 yields the following system of transport equations for the amplitudes: 



(2.8) d t Aj + kj ■ VA 



if a > 1, 

-iV e s{A)Aj if a = 1, 



where we have used the fact that Atpj = 0, in view of (2.7). For a = 1, the effective 
(nonlinear) potential V c s{A) is given by 

J 

V eS (A) :=K*(Y,\Ae\ 2 ). 

i=i 

We see that for a > 1 no nonlinear effects are present in transport equations for 
the leading order amplitudes. The case a = 1 is therefore seen to be critical as far 
as semi-classical asymptotics is concerned. 

Lemma 2.1. The transport equation (2.8) with initial data (oj)j g r £ L 2 A(M. d ) 
admits global-in-time solutions A £ C([0, oo); L 2 D ^l(M d )) ; which can be written in 
the form 

(2.9) Ajit.x) = a 3 (x - tkj)e iS ^' x \ 

where Sj = for a > 1 and Sj G C([0, oo) x M. d ) for a = 1 is given by 

rt J 

(2.10) Sj(t,-) = - [K * y2\ai(-+(T-t)k 3 -Tk e )\ 2 )dT. 

Jo v ' 
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In particular we have mass conservation for each individual mode 

(2.11) ||^(t,-)||£»(R*) = ||o i |U» (R d )) Vf G R. 

In contrast to <f>j the phases Sj are only slowly varying, i.e. they oscillate with 
frequencies larger than 1/s. They describe the nonlinear self- modulation of the 
amplitudes but do not show up in quadratic quantities, like the mass density \Aj\ 2 
etc. 



Proof. By multiplying (2.8) with Aj and taking the real part, we see 

(d t + k r V)\A 3 \ 2 = 0, 

which yields \Aj(t,x)\ 2 = \dj(x — tkj)\ 2 and thus (2.9) and (2.11). Finally, inserting 
(2.9) into (2.8) and integrating along characteristics, we obtain the expression (2.10) 
forSj. □ 

Having obtained the characteristic phases 4>j and the leading order amplitudes 
Aj we shall now turn our attention towards the remainder, i.e. 



(2.12) i?« pp ) 



z 2 X 2 + eY R ifa=l, 
z 2 X 2 + e a {Y + Y R ) ifa>l. 



For a = 1 the term Yr appearing within i?(u| pp ) is formally of order O(e). Thus, 
at first glance, Yr seems to be too large to be considered a part of the remainder. 
It will be our main task to show that Yr is indeed sufficiently small as e — > 0+ . To 
this end, we shall rely on the framework of Wiener algebras. 

Remark 2.2. In the case a > 1, there is a some freedom in deriving the system of 
amplitude equations. Indeed, one could regard e Q_1 merely as a prefactor for the 
interaction kernel K and solve in the next step X\ + e a ~ 1 Y = 0. Instead of (2.8), 
we would consequently get 

d t A) + kj ■ VA) = -vf-^V&iA^A* for a > 1, 

which yields 

(2.13) A?(t,x) = aj (x - tk 3 )e iea ' ls ^ x \ 

with Sj given as before. In other words, the leading order amplitudes would become 
e-dependent and instead of (2.12) the remainder would have the form 

i?« pp ) = e 2 X 2 + e a Y R for a > 1. 

We expect that the subsequent analysis can be adapted to the case of e-dependent 
amplitudes given by (2.13), and in fact we expect a slightly better error estimate 
in this case. Namely, we expect (1.5) to hold with (3 = 1 for all a ^ 1. However, 
we hesitate to carry through this approach in full detail, since the separation of 
slow- and fast scales within the leading order approximation becomes much less 
apparent. In addition, one usually does not regard (2.13) as a proper leading order 
amplitude in semi-classical analysis (due to its e-dependence). We shall therefore 
resume our analysis with e-independent Aj as in our original ansatz (2.1). 
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3. The Wiener algebra framework 

We now present the analytical framework of Wiener algebras which already 
proved its use in similar circumstances, cf. [5, 8, 12]. We start with the follow- 
ing definition. 

Definition 3.1 (Wiener Algebra). We define 

W(R d ) := {/ e S'(R d ; C), ||/|| w := \\f\\ L i m < ^} . 

The following properties of W(K ) have been proved in [8, 12]. 
Lemma 3.2. 

i. W(R d ) is a Banach space, continuously embedded into L°°(R d ). 
ii. W(M ) is an algebra, in the sense that the mapping (f,g) i— ► fg is contin- 
uous from W(R d ) 2 to W(M ), and moreover 

Vf,g€W(R d ), \\fg\\w<\\f\\w\\g\\w- 

iii. Let U e (t) = e l£ i A denote the free Schrddinger group. Then, for all teM, 
U 6 (t) is unitary on W(R d ). 

Assertion iii. follows from the fact that U £ (t) = e 1 *'^ I 2 , acting as a multiplica- 
tion operator in Fourier space. From now on, we shall consider the Cauchy problem 
(1.3) to be posed in W(R d ). To this end we need the following wcll-poscdness result 
(which is an adaptation of the one given in [5] to the case of Hartree nonlincaritics). 

Lemma 3.3. Consider the initial value problem 

(1.3) ied t u £ + —Au E = s a (K * \u e \ 2 )u E , u e \ t=Q = u e , 

with a > 1. If K G L ca (R d ) and u% G W(R d ) with ||«§||w < A), then there 
exists a Tq > 0, which depends on Dq but not on e, and a unique solution u e G 
C([0,T o ];W(R d )) of (1.3). 

Proof. Duhamel's formulation of (1.3) reads 

u £ (t) = U E (t)u Q - ie a - x f U E (t - t) ((A * \u s \ 2 )u e )) (r) dr. 
Jo 

Denote, for fixed ug, the right hand side of this formula by $ e (u e )(i). From 
Lemma 3.2 iii. wc have 

(3.1) \\^(u E )(t)\\ w < Do+e*- 1 f || ((A* K| 2 K)) (r)\\ w dr. 

Jo 

In order to control the nonlinear term, we need to estimate expressions of the form 
(A * (uv))w in W(R d ). To this end, we first use Lemma 3.2 ii, to obtain 

\\(K* (uv))w\\ w ||A* (™)|| w ||w|| w . 

By Holder's inequality we also get 

(3.2) || A'* (uv)\\ w = \\K7^v)\\ L i \\K\\ L ~\\uv\\ w , 
and applying again Lemma 3.2 ii, we arrive at 

(3.3) \\(K * (uv))w\\ w < HAll^lMMMklMk foru,v,weW(R d ). 
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Thus, from (3.1) we obtain 

||$ £ (^)(i)|| w < J D +^- 1 ||^|k» f\\u%T)f w dr. 

Jo 

Moreover, u e i-> <$> £ {u £ ) is locally Lipschitz in U := C([0, T], iy(R d )): If \\u E \\u ^ D, 
\\v £ || u ^ -Dj then there exists C = C(D) such that 

\\^(u £ )(t)-^(v £ )(t)\\ w ^C(D) /Vto-t^Mll^dT, Wg[0,T]. 



A fixed point argument in {h e [/ : sup te j 0T ] ||u(4)||vy ^ I?}, with Z? > Dq, for 
T = Tq sufficiently small, then yields Lemma 3.3. □ 

Having set up an existence result for the exact solution u e in W(R d ), we now 
turn to the approximate solution u| pp given by (1.6). To this end, we shall need 
the following lemma, which shows that the Wiener space is perfectly adapted to 
our problem. 

Lemma 3.4. Let kj g M. d , c, g R, and b g ^(T, W(R d )), and set 

J 

f(x,y) = ^2b j (x)e i ^-y+^. 

3 = 1 

Then, for all e > the function /(•, ■ / e) : x i— » /(x, x/e) /ies m Vy(R d ) wif/i 

||/(-,-A)lk<IHU = EH 6 ill^ 

Proof. We write 

./ J J 

ii/o, -A)ik = II E e< S(- - *i/ e )ILi < E ii^o - ^/ £ )H^ = E IIM^ 1 - 

3=1 3=1 3=1 

The last term is, by definition, ||6||„4,. □ 

Remark 3.5. This lemma in general does not hold for functions of the form 
f(x,y) = J2j=i^ ) j( x ) elipjt ' y ^ with <fj(y) ^ kj ■ y + Cj. A generalization of our 
study to non-plane wave WKB phases therefore seems to be a delicate issue (at 
least within the Wiener framework) and by no means straightforward. 



Since the phases <fij considered in this work are of plane- wave form (2.7), applying 
Lemma 3.4 with Cj = — ^\kj\ 2 to (1.6), we immediately obtain 

(3-4) IK PP («,-)lk<||A(t > -)IU. 

Similarly, we can estimate the expression (2.3) for X2 by 

1 1 3 

(3.5) \\X 2 (t,-)\\w < 2 I|AA(< '' )IU = 2^ AAj{tr)llw - 

3=1 

In addition, we obtain the following estimate for Y, defined in (2.4): 

J 

\\Y(t,-)\\ w ^ \\K*J2\Mt,-)\ 2 \\ \\A(t,-)\\ Al 
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where we have used Lemma 3.4 and Lemma 3.2 ii. This can be estimated further 
similarly to (3.2) by using Lemma 3.4 (with kj = 0, Cj = 0), as well as Assertion 
ii. of Lemma 3.2 and the fact that £°° Cf 1 , to obtain 

(3-6) \\Y(t,.)\\ w ^\\K\\ L ~\\A(t,.)\\ 3 A . 

In order to close the argument, we consequently require appropriate bounds in 
*4(R d ) on the amplitudes Aj(t, •), together with their spatial derivatives. 

Lemma 3.6. Let a ^ 1 and K G L°°(R d ). For all a = (cij)jer G A(R d ), there 
exists a unique solution A G C([0, oo); ,A(R d )) to the system (2.8). Moreover, if 
, i G A(R d ), for \p\ < 2, then (^Aj) j& r G C([0, oo); A(R d )). 

Proof. For a > 1 the statements of the lemma are immediately clear, since in this 
case Aj(t, x) = aj(x — tkj), cf. (2.8). For a = 1 we rewrite (2.8) in its integral form 

(3.7) Aj(t,x)=aj(x-tkj)+ [ N{A)j{t,x + (r - t)kj)dr, 

Jo 

where the nonlinearity J\f(A)j is given by 

a/^), = -*(^*(Ei^i 2 )H- 

Invoking the same arguments as for the derivation of (3.6), we obtain 
\\M{A)\\ A = j2\\N{A)^\ w ^\\K\\ L ~\\A\\\. 

3=1 

This shows that N(A) defines a continuous mapping from .A 3 to A and a local- 
in-timc existence result immediately follows from the standard Cauchy-Lipschitz 
theorem for ordinary differential equations. That the solutions Aj indeed exist for 
all t ^ then follows from the explicit representation (2.9). From the latter we 
additionally obtain the propagation of regularity, by explicit calculation of d§A. □ 

Lemma 3.6 consequently establishes the estimates in VF(M d ) for u| pp , X2 and 
Y in a rigorous way. Note however, that the above given estimates do not yield 
an estimate for the remainder R(u% pp ), given by (2.12), since it also includes Yr, 
which we completely ignored so far. We will make up for it in the following section. 

4. Estimates on the remainder and proof of the main theorem 

It remains to estimate in W(R") the term Yr given in (2.5). To this end we shall 
prove the following key technical result. 

Proposition 4.1. Let Yr be defined by (2.5) with plane-wave phases given by 
(2.7) and assume \A\^f := inf{|/cf — k m \ : £,m G T, i ^ m} > 0. Moreover, let K 
be such that K G L°°(R d ) and VK G L°°(R d ). Then we have the following bound: 

(4.1) \\Y R (t, < eC K \\A(t, OlllflM*. -)IU + l|VA(t, -)|U), 

where Ck > is independent of e. 
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Proof. Recalling the definition of Yr given in (2.5) and taking into account the 
particular plane-wave form (2.7) of the phases tfij, we obtain from Lemma 3.4 and 
Assertion ii. of Lemma 3.2 that 

J 



i(4>i-4> m )/e'\ 



(4.2) \\Y R \\ W < |L4|U|| Y, K*(AiA % 

Using, 

e iyk/e+z = -ie-rjL ■ V y e ty - k / £+z (zeC), 
we can perform a partial integration w.r.t y, and rewrite 

K * (A f I ra e*«-^»/ E ) = [ K(x - y)Mv)A m (.V> ilMy) ~* m(y))/e d V 

V y (K(x - y)Mv)A m (v))e i{Mv) ~* m( - v))/e d V- 



IE 



d \kt - k m \ 2 



To show that the boundary terms vanish, assume first that Aj G S{R d ), the 
set of Schwartz functions (for which the boundary terms clearly vanish). Since 
5(M ) is dense in L 1 (M d ), Fourier transformation implies that 5(M ) is also dense 
in VF(M d ). Consequently, the fact that the expressions on both sides are norm- 
bounded scsquilinear forms establishes the above formula. 

Using kg 7^ k rn G R d , we know that A^ iTO := n^ffrp G R d is well defined and 
consequently 

K * (AtA^'-M/ 5 ) = - ie{k l>m -VK) * (A e A m e l ^-^/ e ) 

+ ieK * ((A f|m -V(4A))e i( *'-* ra)/E ). 
Using the estimate (3.2) and Lemma 3.4 (with J = 1) we get 

\\K*(A e A m e l ^-^/ £ )\\ w s? e\\A itm -VK\\ L oo\\A ( \\ w \\A m \\ w 

+ s\\K\\ L °o || A^ m -V(A e A m )\\ w . 
Invoking again Lemma 3.4 (with kj = 0, Cj = 0) and Lemma 3.2 ii, we conclude 

| K*(A e A m e^-^) 

e,m=i 

.7 

s^e E (d|A^ ro |||Vif|| i oo||^|| W r||yl m || W r + 2||^|| i »|A^ TO |||V^|UP 



rn 1 1 W 



( .111 = 1 



<e|A| 00 p|U(d||Vlir|| i o ||A|U + 2||A-|| L «||VA|| 
where ||V-KT||loo = max ||9 Xrl ii'||x,oo, and 

J J d 

\\va\u = e iiv^iu = EE ll^lk- 

f=l £=1 n=l 

This, together with (4.2) yields the estimate (4.1). □ 
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Remark 4.2. Proposition 4.1 shows that ||Yr||w = 0(e) and thus can indeed be 
considered a part of the remainder. Note that in the proof it is essential to invoke 
a stationary-phase type argument first, before starting to take estimates. In fact 
we would not succeed to show that ||Yr||w = 0(e), if we would estimate Yr in its 
original form. 

By combining the results of Lemma 3.6 and Proposition 4.1, we obtain the 
following result. 

Corollary 4.3. Under the assumptions of Lemma 3.6 and Proposition 4-1, there 
exists, for every T > 0, a constant Cr(T) > 0, independent of s > 0, such that the 
remainder R(u 6 app ) given by (2.12) satisfies 

(4.3) \\R« pp )\\w <^C R Vie[0,r], Ve>0, 

where 

'l !/«e{l}U[2,oo), 
a if a e (1,2). 
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Proof. Lemma 3.6 guarantees the existence of the norms \\A(t, 011-4) II V^4(t, Oll-Aj 
||AA(i, < oo for all t £ [0, oo), which are independent of e > and continuous 
in t. Hence, taking their maximum over t G [0, T], we obtain (4.3) from (3.5), (3.6), 
(4.1) and the definition (2.12) of fi« pp ). □ 

With the estimate of Corollary 4.3 on R(u% pp ) at hand, we can finally state the 
proof of our main theorem, which follows the basic ideas established in [13] for jus- 
tifying the nonlinear Schrbdinger equation as a modulation equation for dispersive 
waves. 

Proof of Theorem 1.2. We consider a fixed T > and introduce the following scaled 
error r e between the original solution u £ to (1.3) subject to the initial data (1.2) 
and the approximation (1.6): 

S r . u — u app , 

with a parameter (3 > to be specified below. Hence, r e (0) = 0. From (3.4) and 
Lemma 3.6 we know that there exists a constant C 'a > 0, independent of e, such 
that |K pp (i, -)\\ w < C A , for all t € [0,T]. Since r £ (0) = 0, it follows \\u e \\ w < C A . 
Consequently, for any D > Ca, Lemma 3.3 yields the existence of a unique solution 
u e £ C{[0,T o ],W{R d )) for some T > with |K(i)lk < D for t e [0,T ]. 
Moreover, from (1.3) and (2.2) it follows that r £ satisfies 

ied t r £ + ^Ar e = e a ^(M(ui pp + £ V) - A^« PP )) - e^R(ut pp ) 

with M(u) = (K * \u\ 2 )u for t ^ r := min{T , T}, and by Duhamcl's formula and 
Lemma 3.2 iii. we obtain 

\\r £ (t)\\ w < gOt—0—i f \\M(ut pp (r)+s^(T))-M(ut pp (r))\\ w dr 

(4.4) Jo t 

+ S-V- 1 I \\R(ul pp (r))\\ w dr, 
Jo 

for alH ^ t. Writing 

M(u + r) — M(u) = (K * (uf + ur + |r| 2 ))(u + r) + (K * \u\ 2 )r, 



12 



J. GIANNOULIS, A. MIELKE, AND C. SPARBER 



the estimate (3.3) gives 

(4.5) ||«M(« + r) - M{u)\\w < ||^||l-(3||u||^ + 3\\u\\ w \\r\\ w + \\r\\ 2 w ) \\r\\ w . 
Hence, replacing u — w| pp and r = e^r s , and recalling j3 > 0, we obtain for any 
C > and s a G (0, 1], such that 3e%C A C + efC 2 = C\, that 

(4.6) \\M(ut pp + eV) - M(ul pp )\\ w ^ e?C m \Y\\ w Ve^o.Kt 

where Cm '■= M\^\\l°°0%, as long as ||r E ||w s% C. 

Inserting the bounds (4.6) and (4.3) into (4.4), and recalling that eq < 1, a ^ 1, 
t ^ T, we consequently obtain for (3 G (0,7—1] 

\\r s (t)\\w ^C r T + Cm I V(T)\\ w d T Ve^e , t^r. 
Jo 

By Gronwall's lemma this yields 

\\r E (t)\\ w ^C R Te c ^ Ve^o.t^r. 

Hence, setting above C := CnTe CMT and D := C + Ca this estimate guarantees 
that the solution u E exists on the whole time interval [0, T], cf. the proof of Lemma 
3.3. Moreover, recalling Lemma 3.2 i, we finally obtain the error estimate (1.5) 
which finishes the proof. □ 

Acknowledgement. The authors want to thank Robert Schubert for pointing 
out to them the alternative approach described in Remark 2.2. 
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